Astronomy <fe Astrophysics manuscript no. Is50-rev2 


© ESO 2008 


February 5, 2008 





GeV-TeV gamma-ray light curves expected in the IC pair 
cascade model for massive binaries: Application to LS 5039 

W. Bednarek 



Department of Experimental Physics, University of Lodz, 90-236 Lodz, ul. Pomorska 149/153, Poland 



o 
o 

(N 
> 

o 



(N 



o 
o 



Received 



accepted 



X 



ABSTRACT 



Context. TeV gamma-ray emission from two massive binaries of the microquasar type, LS 5039 and LS I +61° 303, 
show clear variability with their orbital periods. 

Aims. Our purpose is to calculate the GeV and TeV 7-ray light curves from the massive binary LS 5039 which are 
expected in the specific Inverse Compton pair cascade model. This model successfully predicted the basic features 
of the high energy 7-ray emission from LS 5039 and LS I -|-61° 303. 

Methods. In the calculations we apply the Monte Carlo code which follows the IC pair cascade in the anisotropic 
radiation of the massive star. 

Results. The 7-ray light curves and spectra are obtained for different parameters of the acceleration scenario and the 
inclination angles of the binary system. It is found that the GeV and TeV 7-ray light curves should be anti-correlated. 
This feature can be tested in the near future by the simultaneous observations of LS 5039 with the AGILE and GLAST 
telescopes in GeV energies and the Cherenkov telescopes in the TeV energies. Considered model also predicts a broad 
maximum in the TeV 7-ray light curve between the phases ~ 0.4 — 0.8 consistently with the observations of LS 5039 by 
the HESS telescopes. Moreover, we predict additional dip in the TeV light curve for large inclination angles ~ 60°. This 
feature could serve as a diagnostic for independent measuring of the inclination angle of this binary system indicating 
also on the presence of a neutron star in LS 5039. 
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1. Introduction 

Recently, two massive binary systems of the microquasar 
type have been discovered as TeV 7-ray sources, LS 5039 
(Aharonian et al. 2005) and LS I +61° 303 (Albert et 
al. 2006)). The TeV 7-ray emission is modulated with 
the orbital periods of these two binary systems (Albert 
et al. 2006, Aharonian et al. 2006a). This feature strongly 
indicate that high energy processes occur already inside 
the binary systems not far away from the massive star. 
In fact, such modulation of the 7-ray signals from these 
two binaries (the phases of the maximum emission) have 
been recently predicted based on the calculations of the 
propagation of high energy 7-rays inside massive binaries 
(Bednarek 2006a) . The optical depths for TeV 7-rays inside 
these specific binary systems have been also calculated in 
other papers (e.g. Botcher & Dermer 2005, Dubus 2006). 

In this paper we concentrate on the application of the 
detailed Inverse Compton (IC) pair cascade model for 
the binary system LS 5039 (applied aheady to LS I -1-61° 
303, Bednarek (2006b; B06)). It is usually expected that 
such TeV 7-ray emission can appear in LS 5039 as a result 
of comptonization of the stellar radiation by relativistic lep- 
tons in the jet (e.g. Paredes, Bosch-Ramon & Romero 2006, 
Dermer & Bottcher 2006). Although hadronic production 
of 7-rays has been also considered (e.g. Aharonian et 
al. 2006b, Romero, Christiansen & Orellana 2005). In our 
model, we take into account the effects of anisotropic cas- 
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cade initiated by leptons in the radiation of the massive 
star including the synchrotron losses. We calculate the 7- 
ray light curves and spectra in the GeV and TeV energy 
ranges and discuss them in the context of recently observed 
7-ray light curve from LS 5039 (Aharonian et al. 2006a). 

The basic parameters of the binary system LS 5039 
have been recently reported by Casares et al. (2005): the 
semi-major axis, a = SAr^,., ellipticity e = 0.35 ± 0.04, 
the inclination of the binary system toward the observer, 
i = 24.9 ± 2.8° for the case of the black hole with the mass 
3.7 M0 and ~60° for the neutron star, the azimuthal an- 
gle of the observer in respect to the periastron passage, 
Lu — 225°, radius of the massive star, r^, — 9.3^q'q Kq, and 
its surface temperature, Tg — 3.9 x 10^ K. For these param- 
eters the distance of the compact object from the massive 
star chances in the range: Tp = 2.2rj, at the periastron up 
to Ta = 4.5rj, at the apastron. 

The binary system LS 5039 shows relativistic radio jets 
on miliarcsecond scales, with the speed of ^ 0.3c (Paredes 
et al. 2000). It has been identified to be a counterpart of 
the EGRET source 3EG J1824-1514 which has relatively 
flat spectrum above 100 McV with the spectral index < 2 
(Paredes et al. 2005). Moreover, the position of LS 5039 is 
consistent at the 3(t level with recently detected TeV source 
HESS J1826-148 (Aharonian et al. 2005). The spectrum 
of this source above 250 GeV is also flat with the photon 
index 2.12 ± 0.15, although the flux is about two orders of 
magnitude lower than at GeV energies. Recent re-analysis 
of the TeV 7-ray light curve by Casares et al. (2005), using 
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Fig. 1. The optical depths for 7-rays (as a function of their 
energy) on the pair production in collisions with stellar 
photons. 7-rays arc injected in the jet at the distance from 
its base z = Ir^, ((a) and (c)), and lOr^, ((b) and (d)) for 
the periastron passage (2.2rj,, left figures) and the apastron 
passage (4.5r^, right figures) of the compact object on its 
orbit around the massive star in LS 5039. Specific curves 
show the optical depths up to the infinity (accept the case in 
which photons propagate only to the stellar surface) for the 
injection angles of 7-rays, a, measured from the direction 
defined by the centers of the stars, a = 0° (thick full curve, 
direction toward the massive star), 30° (thick dashed), 60° 
(thick dotted), 90° (thick dot-dashed) , 120° (thin full), 150° 
(thin dashed), 180° (thin dotted). 



new orbital parameters, has shown possible flux variations 
of a factor ^ 3 with the maximum around the phase ~ 0.9. 
In fact, the HESS Collaboration (Aharonian et al. 2006a) 
recently detected 7-ray emission from this source clearly 
modulated with the orbital period of the binary system 
with possibly two peaks at the phases ~0.5-0.6 and ~0.8- 
0.9. 



2. A microquasar in the compact massive binary 

Details of the IC pair cascade model, which we want 
to apply to LS 5039, have been recently described in 
Bednarek (B06, applied for another binary system LS I 
-1-61° 303). In this letter, we remind only the most im- 
portant features of the model, concentrating on the issues 
which might be characteristic for LS 5039. 

Leptons injected into the radiation field of the mas- 
sive star initiate the IC pair cascade provided that, at 
least in some directions, the optical depths for 7-rays (in 
7 -|- 7 ^ e^e~ process) are above unity. We calculate such 
optical depths in the general case, i.e. for an arbitrary place 
of injection of 7-ray photons with arbitrary energies and 
angles of propagation, applying the parameters of the mas- 
sive star in LS 5039 (see Fig.[T]). For certain range of photon 
energies and their injection angles, the optical depths are 
clearly above unity. Thus, the basic conditions for the ap- 
plication of such a model are fulfilled since the electrons in 
the jet preferentially produce the first generation 7-rays in 



directions of the largest optical depth. We consider a sim- 
ple geometrical scenario in which two sided jet is launched 
along the disk axis. For simplicity, we assume that the sur- 
face of the disk is in the plane of the binary system, i.e. 
the jet direction is perpendicular to the plane of the bi- 
nary system. Relativistic leptons are injected along the jet 
at some range of distances from its base, starting at Zmin 
up to Zmax- The magnetic field strength, which determine 
the process of acceleration of leptons in the jet and their 
losses on synchrotron process, is described by a parameter 
T] (the ratio of the magnetic field energy density to radia- 
tion energy density in the inner part of the accretion disk, 
see details of the model for the magnetic field in the jet in 
Bednarek (B06)). The energy losses of leptons on IC scat- 
tering in the Thomson and the Klein-Nishina regime are im- 
portant as well. They are taken into account when consid- 
ering the acceleration of electrons. We assume that leptons 
are accelerated with a power law spectrum (and spectral 
index which does not depend on z) in multiple shocks prop- 
agating along the jet. The maximum energies of electrons 
are determined by the acceleration mechanism (defined by 
the acceleration parameter ^) and the energy losses on the 
most efficient radiation process (see B06). These maximum 
energies are shown for specific values of ^ and 77 as a func- 
tion of the distance, z, from the base of the jet (Fig.[2|). It is 
clear that for reasonable values, leptons can reach energies 
over a few tens of TeV, already at distances from the base of 
the jet closer than 10 radii of the massive star. We consider 
two different scenarios for the injection rate of relativistic 
leptons along the jet. The injection rate (electrons per unit 
length along the jet per unit time) is constant along the jet 
from its base at Zmin up to ^imax (i-C N[z) oc const), or it 
drops along the jet proportionally to the magnetic energy 
density in the jet starting from some distance, Zmin, from 
the base of the jet (i.e. N{z) cx z~^). 

3. The cascade 7-ray spectra 

Since leptons are injected isotropically in the jet, they pro- 
duce 7-rays preferentially in directions where they find the 
largest optical depth. Therefore, absorption of 7-rays and 
interaction of secondary pairs have to be taken into ac- 
count when calculating the 7-ray spectra produced by lep- 
tons in the radiation field of the massive star both at the 
periastron and the apastron passages. We have performed 
the Monte Carlo calculations of the 7-ray spectra escaping 
from the binary system at arbitrary directions, applying 
the above mentioned model for the acceleration of primary 
leptons in the jet. To allow a simple analysis, at first we 
consider the case of injection of leptons in the jet with the 
power law spectrum, and differential spectral index equal 
to —2, with the cut-off which depends on the location of 
the acceleration place in the jet and on the phase of the 
compact object (see Fig. 2). These cut-offs are determined 
by the values of 77 and Since two jets are expected from a 
single compact object, we consider the production of 7-rays 
in the jet which propagate above the plain of the binary sys- 
tem (i.e. the jet directed towards the hemisphere containing 
the observer) and in the jet propagating below the plane of 
the binary system (i.e. the counter-jet), see (B06). The re- 
sults of example calculations of the 7-ray spectra from the 
IC pair cascade are shown in Fig.[3]for different values of 
^ and 77. The observer is located at the inclination angles of 
25° and 60°. The spectra are shown for the compact object 
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Fig. 2. The maximum Lorentz factors of electrons as a func- 
tion of the distance, z, from the base of the jet (measured 
in units of the stellar radius) for different acceleration ef- 
ficiencies: ^ = 0.5 (from the upper, thickest curve), 0.1, 
0.03, and 0.001 (to the bottom, thinner curve) and the pa- 
rameter describing the magnetic field at the base of the 
jet (at at the inner radius of the accretion disk) equal to 
B = 2> X 10^ G (a) (77 = 1 - the equipartition of magnetic 
field with disk radiation), 77 = 0.1 ((b) and (d)), and 0.01 
((c) and (e)). The Maximum Lorentz factors are obtained 
from the comparison of the acceleration rate with the en- 
ergy loss rate on synchrotron process (dotted lines), inverse 
Compton process in the Thomson regime (full curve) and 
the Klein- Nishina regime (dashed curves), for the pcrias- 
tron (figures (a), (b), and (c)) and the apastron ((a), (d), 
and (e)) passages of the compact object. 

at the periastron (left figures) and the apastron passages 
(right). The 7-ray spectra, produced by Icptons injected 
in the jet at the periastron passage, steepen significantly 
at energies above a few tens of GeV due to the efficient 
cascading. This cascading effects are significantly smaller 
when the compact object is close to the apastron. However, 
note that the 7-ray spectra calculated for these parame- 
ters do not fit nicely to the spectra observed recently by 
Aharonian et al. (2006a) due to the deficit of 7-ray emis- 
sion above ~ 100 GeV. Therefore, in the next section we 
consider the model with the fiatter profile for the injection 
of electrons into the jet (i.e. N{z) = const). 

4. Confrontation with observations of LS 5039 

Let us select the set of parameters which allow acceleration 
of electrons to energies up to '^lO TcV (see upper panels 
in Fig. [2]). For these parameters (^ = 0.1 and 77 = 0.1) 
and the case of electrons injected with the constant rate 
along the jet between Zmin = O.lr^ and Zmax = 10^,^, we 
calculate the 7-ray light curves in two energy ranges, i.e 
1-10 GeV (the GeV region, the AGILE and GLAST en- 
ergy range) and >300 GeV (the TcV region, the Chcrcnkov 
telescopes energy range), and two possible inclination an- 
gles of the binary system LS 5039, i.e. 25° and 60° (the 
range of angles derived by Casares et al. 2005). The TeV 
7-ray light curves, calculated for 25°, show broad maximum 
between the phases ^ 0.4 — 0.8 peaked at the phase ~ 0.6. 
On the other hand, the light curves for the inclination angle 
60° show more complicated structure. The broad maximum 
between the phases ^ 0.4 — 0.8 shows clear two additional 
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Fig. 3. Spectral Energy Distributions (SED) from cascades 
initiated by primary electrons injected in the jet with: the 
efficiency of electron injection depending on the distance 
from the base of the jet as N{z) oc (independently on 
the phase of the binary system), and differential power law 
spectrum of electrons with the index equal to —2. The left 
figures show the 7-ray spectra, produced at the periastron 
passage of the compact object for the observer located at 
the inclination angles i = 25° and 60°, from the jet and 
counter-jet (the sum of both) . The jets propagate perpen- 
dicular to the plane of the binary system. The 7-ray spectra 
produced at the apastron passage of the compact object are 
shown on the right figures. The specific spectra are calcu- 
lated for the acceleration conditions in the jet described by 
the parameters: ^ = 0.03 and 77 = 0.1 (full curves), ^ = 0.3 
and 77 = 0.01, (dashed), and ^ = 0.3 and 77 = 0.1 (dot- 
dashed) . 
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Fig. 4. The 7-ray light curves in the energy range 1 — 10 
GeV (dashed curves) and >300 GeV (full curves), produced 
by electrons in the jet (thin curves) and the jet -I- counter- 
jet (thick curves), are shown on the upper panel for two 
possible inclination angles of the binary system LS 5039 
(25° and 60°), ^ = 0.1, 77 = 0.1, z^i^ = O.lr,, and z^,^ = 
lOr^,. On the bottom panel, the 7-ray light curves are shown 
for ^ = 0.2, 77 = 0.1, and z,nax = 20r*. 
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peaks centered at the phases ^0.5 and ~0.8 and a dip at 
the phase ~ 0.65 (Fig. [?]). The general broad maximum 
observed for both inclination angles is consistent with re- 
cently reported TeV 7-ray light curve observed by the HESS 
Collaboration (see the bottom panel in Fig. 5 in Aharonian 
et al. 2006a). Note however, that the HESS TeV 7-ray hght 
curve also shows, although not statistically very significant, 
a dip in the maximum at the phase ~ 0.65, i.e. correspond- 
ing to the interior conjunction of the compact object. The 
confirmation of such a structure by the future HESS obser- 
vations will strongly argue for the inclination of the binary 
system closer to 60°, supporting the hypothesis on the pres- 
ence of a neutron star in this binary system. 

Note, that recently Bottcher (2006) argues for the in- 
clination angle of the LS 5039 closer to 30°, based on the 
deabsorption of the 7-ray spectrum observed by the HESS 
experiment from LS 5039. However, these calculations base 
on the assumption that the cascading processes initiated by 
secondary leptons (originated in the absorption process of 
TeV 7-rays) in the radiation of the massive star can be ne- 
glected. This is valid provided that the magnetic field inside 
the volume of the binary system is strong enough to guar- 
antee efficient cooling of secondary pairs on synchrotron 
process (in the case of LS 5039 should be of the order of 
- 100 G). In another CcLS6, clS discussed in this paper, the 
IC pair cascades should be taken into account since the 
cross sections for 7-7 absorption and IC scattering in the 
KN regime are similar, see e.g. the optical depths for these 
two processes in Fig. 2 in Sierpowska & Bednarek (2005). 

In our model the broad maximum between phases ~ 
0.4 — 0.8 in the TeV 7-ray light curve is due to the favorite 
location of the compact object in respect to the massive 
star and the observer. At these phases the jet is generally 
on the other site of the massive star in respect to the ob- 
server and the conditions for scattering of stellar radiation 
by electrons in the jet into the TeV 7-ray energy range are 
optimal. However, when the line of sight to the compact 
object passes very close to the massive star, which occurs 
in the case of large inclination angles, absorption of TeV 
7-rays in the stellar radiation start to play dominant role. 
This is the reason of the appearance of additional dip in 
the TeV light curve at phase ~ 0.65 for large inclination 
angles of the binary system. For low inclination angles, the 
line of sight is sufficiently away from the massive star and 
the absorption effects of TeV 7-rays are too low to produce 
additional dip in the light curve. 

On the other hand, the GeV 7-ray light curve behaves 
very different. It is clearly anti-correlated with the TeV 7- 
ray emission. This is caused by more efficient cascading 
close to the periastron passage of the compact object in 
which many GeV photons are created. In directions where 
the optical depths of TeV 7-rays are large, the energy from 
electrons, originally transfered to primary TeV 7-rays, is 
degraded to the GeV 7-rays. In contrary, in directions 
of low optical depths, TeV 7-rays escape with relatively 
mild absorption producing the peak in the TeV light curve. 
Therefore, the maximum in the GeV 7-ray light curve ap- 
pears just after periastron, at the phase ~0.05 (i.e. at the 
superior conjunction of the compact object), and the min- 
imum at the phase ~0.7 (i.e. at the inferior conjunction). 
It is expected that the GeV 7-ray emission may vary by a 
factor up to ~30 with the period of LS 5039 (see Fig. [4]). 
Hence, we predict also strong modulation of the GeV emis- 



sion from LS 5039 with its orbital period which should be 
easily observed by the AGILE and GLAST detectors. 

Note that the TeV 7-ray light curve, reported by 
Aharonian et al. (2006a, Fig 5), show clear base line emis- 
sion also at the periastron passage of the compact object. 
In order to check if this emission can originate farther along 
the jet, we show in the bottom panel in Fig.|3]the 7-ray light 
curves for the case of injection of electrons in the jet up to 
^max = 20r*, ^ = 0.2, and t] = 0.1. In fact, such base line 
emission at a level of ~20% of the peak emission appears 
close to the periastron for the inclination angle i = 60°. 
Moreover, previously found two peak structure in the TeV 
7-ray light curve (or broad high level emission through the 
phases ~0.4-0.8) for the case of inclination angle of the bi- 
nary system i = 60° is still present. However, in contrary to 
the previously discussed case in the upper panels in Fig. 4, 
the TeV 7-ray light curve for the inclination angle 25° is 
quite flat due to similar conditions farther along the jet in- 
dependently on the periastron or apastron location of the 
compact object. 

The spectral features observed by the HESS (Aharonian 
et al. 2006a), flat power spectrum at the inferior conjunc- 
tion and steep spectrum at the superior conjunction, can 
be also explained in terms of the considered here general 
model. Such features appear even in so simplified model for 
the injection of rclativistic electrons in the jet (see Fig. 3), 
i.e single power law spectra for the electrons and their injec- 
tion rate independent on the phase of the compact object. 
Based on these results we conclude that a more realistic 
model (which relax these simple assumptions) can certainly 
produce exact fitting of the experimental spectral data, al- 
though huge computation time requirements prevent such 
direct comparison at present. Since in the present model 
the cooling of electrons in the jet on IC process is very effi- 
cient (Bednarek 2006b), and the injection rate of electrons 
does not depend on the phase of the compact object, the 
total requirement on the energy budget of accelerated elec- 
trons is comparable to the power emitted in the 7-rays at 
the high state, i.e. close to the inferior conjunction. 

Acknowledgements. I would like to thank the anonymous referee for 
valuable comments. This work is supported by the Polish MNil grant 
No. 1P03D01028. 



References 

Aharonian, F. et al. (HESS coUab.) 2005 Sci. 309, 746 

Aharonian, F. et al. (HESS coUab.) 2006a A&A, in press 

( astro- ph/0607192) 
Aharonian, F., Anchordoqui, L., Khangulyan, D., Montaruli, T. 2006b 

J.Phys.Conf.Ser. 39, 408 
Albert, J. et al. (MAGIC coUab.) 2006 Sci. 312, 1771 
Bednarek, W. 2006a, MNRAS 368, 579 

Bednarek, W. 2006b, MNRAS 371, 1737 ( B06) 

Bottcher, M. 2006, APh, submitted ( astro-ph/ 0609136} 

Bottcher, M., Dermer, CD. 2005, ApJ 634, L8T^ 

Casares, J., Ribo, M., Ribas, I., Paredes, J.M., Marti, J., Herrero, A. 

2005 MNRAS 364, 899 
Dermer, CD., Bottcher, M. 2006, ApJ 643, 1081 
Dubus, G. 2006 A&A 451, 9 

Paredes, J.M., Bosch-Ramon, V., Romero, G.E. 2006, A&A 451, 259 
Paredes, J.M., Marti, J., Ribo, M., Massi, M., 2000 Sci. 288, 2340 
Romero, G.E., Christiansen, H.R., Orellana, M. 2005 ApJ 632, 1093 
Sierpowska, A., Bednarek, W. 2005, MNRAS 356, 711 



